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ABSTRACT 



Oh' 



Context. The use of hydrodynamical simulations, the selection of atomic data, and the computation of deviations from local thermo- 
dynamical equilibrium for the analysis of the solar spectra have implied a downward revision of the solar metallicity. We are in the 
process of using the latest simulations computed with the C05BOLD code to reassess the solar chemical composition. Our previous 
Ql, analyses of the key elements oxygen and nitrogen have not confirmed any extreme downward revision of Z. 

, Aims. We determine the solar photospheric carbon abundance by using a radiation-hydrodynamical C05BOLD model, and compute 

the departures from local thermodynamical equilibrium by using the Kiel code. 

Methods. We measure equivalent widths of atomic C i lines on high resolution, high signal-to-noise ratio solar atlases of disc-centre 
intensity and integrated disc flux. These equivalent widths are analysed with the use of our latest solar 3D hydrodynamical simulation 
computed with C05BOLD. Deviations from local thermodynamic equilibrium are computed in ID with the Kiel code, using the 
5^ ■ average temperature structure of the hydrodynamical simulation as a background model. 

Results. Our recommended value for the solar carbon abundance, relies on 98 independent measurements of observed lines and is 
A(C)=8.50±0.06, the quoted error is the sum of statistical and systematic error. Combined with our recent results for the solar oxygen 
and nitrogen abundances this implies a solar metallicity of Z = 0.0154 and Z/X = 0.021 1. 

Conclusions. Our analysis implies a solar carbon abundance which is about 0.1 dex higher than what was found in previous analysis 
based on different 3D hydrodynamical computations. The difference is partly driven by our equivalent width measurements (we mea- 
sure, on average, larger equivalent widths with respect to the other work based on a 3D model), in part it is likely due to the different 
properties of the hydrodynamical simulations and the spectrum synthesis code. The solar metallicity we obtain from the C05BOLD 
analyses is in slightly better agreement with the constraints of helioseismology than the previous 3D abundance results. 

VO ■ 

^s^j ■ Key words. Sun: abundances - Stars: abundances - Hydrodynamics - Line: formation 

O 

' 1. Introduction ses of the photospheric solar carbon made in 1980-2000, were 

t-H . converging toward the value of A(C)=8.52 +0.06 (Grevesse & 

■ The importance of an accurate knowledge of the solar abun- Sauval. 119981) , which was slightly lowering the previous values 
. ^ | dances can hardly be overstated since they serve as the reference by including the appropriate NLTE corrections. 

K> for all other celestial objects. The high performance of the new- However, lAllende Prieto et alj d2002h announced an impor- 

■ generation instruments allows to derive accurate stellar abun- tant downward revision of the C abundance from the analysis of 
03 dances and therefore the requested accuracy of the reference tne forbidden [ CI1 872.7 nm line (A(C )= 8.39 + 0.04). A subse- 

solar abundances is increased. This can, at least partly, explain que nt paper bv lAsplundet all d2005al) obtained a similar down- 
fhe current revival in spectroscopic solar abundance studies. The war d revision of the carbon abundance, also when using permit- 
very large gap in resolution between solar and stellar spectra, t ed atomic and molecular lines. 

that existed until a few decades ago, is diminishing rapidly. The The ne w abundances of C, as well as those of other elements, 

majority of recent solar abundance determinations relies on ob- ar e in conflict with some solar properties: solar models ( Yang & 

servational data that are almost 30 years old, both for the disc- Bi. 120071) and helioseimology dBasu & AntiaL 12008b Chaplin & 

centre intensity and for the integrated disc flux (Jungfraujoch Basu. 120081: iDelahave & Pinsonneauli 120061) cannot be recon- 

grating spectra and Kitt Peak Fourier Transform Spectra, respec- di ed with the recent revision of sol ar abundances bv Asplund et 

lively). a i. ( 2005b). 

For a long time, solar abundances were considered well es- Solar abundances of the light elements, which have the high- 

tablished, and only minor refinements were suggested by each est cosmic abundance, are particularly important to understand 

new study, usually driven by improved atomic or molecular data, stellar and galactic composition. Besides being the main contrib- 

By using atomic or molecular lines, or both, the many analy- utors to the solar metallicity Z, the CNO abundances are useful 
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to study the depletion in the interstellar gas (ISM). For example, 
the comparison of the C/O ratio in the ISM with that in the so- 
lar photosphere tells us how much C has been locked into dust. 
Also the study of diffusion effects in the corona and solar wind 
requires the use of photospheric solar abundances as a reference. 
Carbon, being a highly volatile element, partly escaped from car- 
bonaceous chondrites, so that the solar system abundance of C 
relies mainly on the analysis of the photospheric spectrum. 

The solar spectrum is rich in atomic C lines as well as in lines 
of C-bearing molecules. Due to its high first ionisation potential 
(11.26eV), the measurable lines of carbon in the Sun are only 
those of C i. Several tens of C i lines are present in the visual 
and near IR spectrum, but only few are suitable for abundance 
analysis. The chosen lines should be weak, unblended, with ac- 
curately known transition probabilities and, ideally, formed in 
LTE. Strong lines, with large equivalent width (EW>15-20pm) 
should be rejected because the collisional damping constants are 
uncertain. Only one forbidden line, at 872.7126 nm, has been 
detected in the solar spectrum. 

Molecular lines are highly temperature sensitive and require 
a very accurate a nalysis of the pho tospheric thermal structure as 
the o ne made by Avres et all d2006h for the infra-red CO fea- 
tures. Holwege rl OOP! " prefers to consider only atomic lines, 
while Grevesse et alj (1 19871) derive the abundances from the 
vibration-rotation and pure rotation lines of the CO and CN di- 
atomic molecules to be more accurate. 

In the present paper we analyse only C i atomic transitions to 
derive the solar photospheric carbon abundance. 



2. Selection of lines 

As a starting point, we looke d at a sample obtained by com- 
bining the C i lines chosen b vlGrevesse et alj d!991h. Bie mont 
et al. d!993h . iTakedal (11994 . and lAsplund et all (l2005ah (see 
TableQ]). We examined these lines, compared the available solar 
atlases among them and to synthetic profiles. We excluded from 
our analysis the lines that we judged too heavily blended (e.g. 
the line at 477.0 nm) compared to the synthetic spectra, or the 
lines for which the disagreement among observed data was un- 
explained and too large (e.g. the line at 1 180. 1 nm). Furthermore, 
we eliminated the lines whenever we suspected a significant con- 
tamination from telluric absorption, based on the comparison of 
the observed atlases and synthetic spectra, and also on the in- 
spection of spectra of fast rotating stars indicating the presence 
of telluric lines (e.g. the line at 1602.1 nm). The excluded lines 
are flagged by "3" in the columns "Quality" of Table Q] The fi- 
nal list of our sample of lines, labelled as Quality "1" or "2" in 
Table Q] is given in Table [3] We labelled as "1" the lines that are 
not blended, or the blends are negligible in comparison to the C i 
line, or we think we are able to model the blends. We labelled 
as "2" the lines we are less confident in. These lines show dif- 
ferences in the observed spectra (e.g. the line at 7 1 1 . 1 nm) or we 
can hardly reproduce their shape with a synthetic profile (e.g. 
the line at 1734.6 nm) or we are not confident to be able to take 
into account the telluric absorptions (e.g. the line at 1778.9 nm). 
The final selection consists of 45 individual lines for which we 
have 98 EW measurements. The subsample of good data, la- 
belled "Quality=l", contains 25 lines, 66 measured EWs. 

For the abundance determination one could rely on line pro- 
file fitting or on EW measurements. The line profile fitting pro- 
cedure has many advantages due to the fact that not only the 
strength of the line is taken into account, but also the line shape. 
In case the synthetic line profile provides a faithful reproduction 
of the line shape, we consider this procedure superior. We stress 



Table 1. Lines considered for the abundance determination. 
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Quality: 1 good line, 2 line with problems, 3 line rejected 



here that by "line profile" fitting we mean fitting with a syn- 
thetic profile, computed using all the known lines in the range. 
But "fitting" with a synthetic profile consisting of a single line 
is conceptually identical to measuring the EW by fitting with a 
Gaussian or Voigt profile, although it has the advantage of treat- 
ing correctly the line asymmetry which, however, is in general 
irrelevant for abundance work. If poorly known blends interfere, 
the EW measurement procedure with deblending (see below) is 
the more secure option. The present analysis is based on EW 
measurements. We give preference to this approach due to the 
following problems with the available C i lines: 

- a large fraction of the lines are blended, and the atomic data 
of the contaminants are not well known, so that when in- 
cluded in the 3D synthetic spectra, the comparison with the 
observed spectra is not reliable; the measurement of the EW, 
on the other hand, can be reliable, since the extra absorption 
can be modelled by a suitable gaussian or Voigt profile; 

- some of the lines are contaminated by telluric absorption; 
also in this case, the contaminating telluric absorption can be 
modelled as above providing a reliable measure of the EW; 

- for some lines the continuum placement is problematic, due 
to the presence of neighbouring lines, whose atomic data are 
often poorly known; the EW measurement with splot, on 
the other hand, is designed to handle such situations; 
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- NLTE effects, not taken into account in the 3D synthetic pro- 
file, can change the shape of the line. 

The last point has not yet been investigated in detail since no 
3D-NLTE analysis for carbon is available at the moment. But 
we expect that carb on does not behave d ifferently from oxygen, 
and it is shown in lAsplund et all (12004) that the 3D-NLTE line 
profile is different from the 3D-LTE one. One could consider to 
use the line profile fitting technique for clean lines, that form 
a small subsample of the complete set of lines, reserving the 
EW measurement to the "problematic" lines. However, in this 
way the analysis would not be homogeneous over the complete 
sample of lines. The adopted method of EW measurements also 
allows a more direct comparison with other analyses available in 

the literature. 

For the measurement of the EW we used the IRAtP HTodvL 



Table 2. Comparison of the disc -centre EW of two C i lines as 
determined by different authors. 



1993) task splot. In the case of blended lines we used the de- 
blending option of splot, that permits to fit the spectral pro- 
file with a number of Gaussian and/or Voigt functions. In this 
way any known line in a range can be simulated with a theoret- 
ical profile. Generally for weak lines we used a Gaussian pro- 
file to fit the observed profile, while for strong lines we used a 
Voigt function. For unblended lines we also used direct integra- 
tion. We are aware that the observed profile is asymmetric while 
both Gaussian and Voigt functions are symmetric. Several ex- 
periments convinced us that the use of a Voigt profile to measure 
the EW of an asymmetric 3D profile differs from the real EW 
by less than 1 %. This error is surely negligible when compared 
to the uncertainty of the EW measurement due to the continuum 
placement, which in the case of a typical observed spectrum can 
be up to 5 %. 

For the majority of the li nes the equivalent wi dths we ob- 
tain are close to the values o f Biem ont et alj d!993h . but not for 
all. We could compare only the Delbouill e disc-centre spectrum , 
which is the observed data considered in IBiemont etal] dl993). 
In principle, strong lines should be rejected because of uncer- 
tain values of NLTE corrections and line broadening parameters 
which become important. We keep these strong lines anyway in 
the sample, because they do not disagree with the other lines, 
and there is no evident trend of the abundance as a function of 
the equivalent width. 

W hen available, we used log gf from NIST (IWiese et al 
119961) . as retrieved from the ASD database dRalchenkd 12005 
The values ar e given in Table [3] All log gf -values used in 
Biem ont et alJ d 19931) . except the one of the 801.8 nm line, are 
very close to the values of NIST. For o ur sample of Ci lines 
the NIST database relies on four sources ( Luo & Pradhan, 1989; 
iHibbert et aUll993tlNussbaumer & StorevLll984tlWeissUl996h . 
the main one being IHibbert et aTT i 1993 ) which covers all the 
lines. 

For the Van derWaals broadening constants we proceeded as 
in lCaffau et al.1 d2008h. When availab le (for 35 lines of our sam- 
ple) we rely on Barklem et alJ dl998h values. For the remaining 
lines we used the WIDTH approximatio n, implemented in the 
Kurucz routine WIDTH (see Rvan, 1998). If we re move the lines 
withou t Van der Waals broadening constants from Barkler rTet alJ 
dl998l) . the derived carbon abundance is less than 0.03 dex higher 
than when considering the complete sample of lines. Therefore 
we decided to keep all the lines we selected for the abundance 
determination. 
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1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



Notes: G91:lGrevesse et alJ dl991h . B93: IBiemont et ail dl993l) . A05: 
Asplu nd et al.1 (12005a). and C10: this work. 



3. Equivalent widths in the literature 

Reliable observed EWs are those measured bv Grevesse et 
al. d 19911) and by IBiemont et al.l d 19931) on the disc-centre 
Jungfraujoch Atlas. The values used bv lStii renburg & Holweger 
(1990) are taken from the EW measurements bv Baschek & 
H olweger d 19671) which w ere based on old atlases, and the ones 
in I Asplund et al.l d2005ah are the EWs of the synthetic best fit 
profile. Only two lines have been considered in all four analy- 
ses, and they are presented in Table|2] 

For the line at 960.3036 nm, we find a significantly larger 
EW than the other three authors, whose results agree closely. On 
the other hand, our EW for the line at 210 2.3151 nm is very sim - 
ilar to those bv lGrevesse et al.1 dl991h and Biemo nFet al.l d!993l) . 
while the theoretical EW derived bv lAsplund et alJ (|2005a) from 
their 3D model by using their best fit abundance is much lower 
(see Table even though the measured EWs should be cor- 
rected for blending; 

We note that Bie mont et al.l ( 1993) gave a low weight to 
both of these lines, presumably because they are affected by tel- 
luric absorption and other blends, and hence reliable equivalent 
widths are difficult to measure. Nevertheless, Table [2] demon- 
strates once again that equivalent width measurements differ 
considerably from author to author and are a major source of 
uncertainty. 

All the investigations of the solar carbon abundance cited 
ab ove rely on a sing le solar atlas. In fact, as already mentioned 
in lCaffau et al.l d2008l) . the available solar atlases do not always 
agree. This could be due to telluric absorption, to variability in 
solar spectrum, or to systematic effects related to the different 
observations. The present analysis is based on four different so- 
lar spectra and hence should yield more reliable abundances. 



4. Observed spectra 

We considered the same four obse rved solar atlases pu blicly 
available that we already used in Caffa u et al.l d2008l) . For 
disc-centre, this is the double-pa s s gra ting spectrum taken at 
Jungfraujoch by Delboui iFeFaT] dl973l) . ranging from 300 to 
1000 nm, and the infrar ed FTS spectrum taken at Kitt Peak by 
Delbouill e~et~aT] d!981l) . covering the wavelength range 1000 
to 5400 nm (together called Delbouille intensity, DI ). In addi- 
tion, th e disc-centre FTS spectrum published by Nec kel & Labsl 
d!984l) is used (330 to 1250 nm. Neckel intensity. NT). Neckel & 
Labs dl984h : lNeckell dl999l) also provide a corresponding FTS 
spectrum for the integrated disc flux (Neckel flux, NF). Another 
set of Kitt Peak F TS scan s by Br ault and Testerman has been 
made available bv iKuruca d2005al) (300 to 1000 nm, hereafter 
Kuru cz flux, KF). For on e line we resor ted to the ATMOS solar 
atlas dFarmer et al.Lll989l:lFarmerLll994 . 
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5. Model atmospheres 

We used the same radiation-hy drodynamical model , com- 
puted with the CO 5 BO LD code dFrevtag et al.L 120021 120031: 
IWedemever et all |2004) used in our previou s solar abundance 
analys e s. Details of the mod el can be found in Caffau & Ludwig 
d2007l) : ICaffau et all d2008l) . The same holds for the employed 
ID models. As a reference model we used a plane parallel 
ID model (IDlhd with mixing-length parameter of 1.0) that 
shares the micro-physics and radiation transfer scheme with 
C0 5 BOLD. We also used a ID model obtained by temporal and 
horizontal averaging of the 3D hydrodynamic structure on sur- 
faces of equal optical depth (( 3D)), as well as the Holweger- 
Miiller semi-empirical mod el dHolweeerl 1 19671: Holweger & 
Miiller!l97ir When necessary this was put on a column mass 
scale, assuming the same abundances and opacities as in the 
C0 5 BOLD model. For the spectrum synthesis based on the ID 
models, we adopted a micro-turbulence of 1.0 km s _1 , both for 
disc-centre (intensity) and integrated disc (flux) spectra. 

6. Results 

6.1. The [CI] line at 872. 7 nm 

There is only one observable forbidden [CI] line, located at 
872.7126nm (2p 2 : D2 - 2p 2 'So), with a lower level excita- 
tion potential of 1.264eV. This line is important, being weak 
and therefore insensitive t o the assumption about the dam ping 
constant, and according to IStiirenburg & Holweger! d!990h not 
affected by NLTE. 

We measured the EW on the two disc-centre and two 
integrated disc solar spectra. The result is EW(DI)=0.51 1 pm, 
EW(NI)=0.509pm, EW(NF)=0.517pm, EW(KF)=0.508 
pm. We subtracted the contribution of the Fei blending 
line (,1=872.71 32 nm, logg/=-3.93, x =4.186eV, EW(Int)= 
0.040 pm, EW(Flux)=0. 045 pm) according to the 3D computa- 
tion. The [CI] line is formed in LTE, so that the LTE abundance 
derived from the four observed solar atlases should be in close 
agreement. While the two disc-centre and the two integrated 
disc spectra are in very good agreement with each other, we 
find a difference in the carbon abundance of 0.04 dex between 
disc-centre and integrated disc spectra. This effect could be 
attributed to NLTE effects on the blending iron line which can 
be different for disc-centre and integrated disc spectra. 

6.2. Permitted lines 

The abundances derived from all the lines of our sample have 
been assembled in Table [3] Atomic data for each line, the mea- 
sured EWs, and the derived LTE abundances using the 3D, (3D), 
and IDlhd models are provided in the table. The total 3D cor- 
rection, defined as 3D - IDlhd is positive for all the lines, except 
for the 477.5907 nm line. The so-called granulation correction, 
quantified by the difference 3D-(3D) measures the effect of the 
horizontal temperature fluctuations. It is never very large, and 
may be in the same direction as the total 3D correction, or in 
the opposite direction, depending on the line. All our permitted 
Ci lines originate from highly excited lower levels (x > 7 eV). 
The 3D-(3D) corrections for the weaker lines (EW < 4.0 pm) in 
the optical and near infra-red range (450 nm< A <820 nm) are 
small and mostly negative (3D-(3D)< 0.01). We can compare 
these 3D-(3D) corrections to the value A„ ra „ defined in Steffen 
& Holweger (2002). The hydro-simulation they consider is a 2D 
model, while the ID model is the temporal and horizontal av- 



erage over surfaces of equal optical depth of their 2D model, 
meaning t hat it is similar to our (3D) model. The C i lines inves- 
tigated in Ste ffen & Holweger] d2002l) have /1=550 nm and EW« 
0.5 pm. Forx > 6eV, A glan is indeed slightly negative, (see their 
Table 1 and their Figure 5) in qualitative agreement with the re- 
sults found for comparable C i lines in the present study. 

Both 3D corrections (3D -IDlhd and 3D-(3D)) increase 
with EW, possibly indicating an inadequate choice of the micro- 
turbulence parameter for the ID models. We cannot discern any 
trend with the excitation energy or the wavelength. 

We do not yet have the capability of computing the devia- 
tions from local thermodynamic equilibrium (NLTE effects) in 
the 3D spectrum synthesis, and, to our knowledge, such calcu- 
lations have not yet been performed elsewhere. As a first ap- 
proximation we have therefore computed ID NLTE corrections 
using the (3D) model as a background model. For each line we 
computed NLTE corrections with the Kiel code (Steenbock & 
Holweger, Il984h and the model-atom of Sturenburg & Holwegerl 
d!990l) . The line blanketi ng is treated with an opac ity distribution 
function as provided bv lCasteili & Kuruca (12003). assuming so- 
lar metallicity and a micro-turbulence of 2kms~' . We consid- 
ered three possible choices for the parameter Sh quantifying the 
thermalizing effect of collisions with neutral hydrogen acco rd- 
ing to the ge neralised Drawin approximatio n dDrawinL 119691) as 
proposed by Steenb ock & Holwegerl d!984l) : 

1. classical scaling (Sh = 1); 

2. no effect of collisions with neutral H (Sh = 0); 

3. intermediate collisional efficiency (Sh = 1/3). 

The NLTE correction obtained for each line is listed in 
Table[3] The re is a general good agreement w ith t he NLTE com- 
putation s of IStiirenburg & Holwegerl (11990) and Asplund et ID 
(I2005al) . with a maximum difference of 0.02 dex. We recall that 
these studies rely on a different model atmosphere for the NLTE 
computation, so that differences of a few hundredth of dex can 
easily be attributed to the different input solar model. This agree- 
ment in the 1D-NLTE computations is encouraging. As long as 
no 3D-NLTE computation is available, it is certainly justified to 
apply this 1 D-NLTE correction to our 3D-LTE abundances, as 
also done by Asplu ncTet alj d2005al) in their careful work. 

As explained above, we rely on EW measurements for the 
carbon abundance determination, because a considerable frac- 
tion of the lines is blended, and NLTE effects are non-negligible. 
Nevertheless, it is useful to compare the observed spectrum with 
the 3D synthetic profiles for some of the cleanest lines. A few ex- 
amples are shown in Fig.Q] The agreement is encouraging, but 
the abundance needed to achieve the best (visual) agreement be- 
tween 3D-synthetic and observed line profile is not always iden- 
tical to the abundance obtained from matching the EWs. This can 
be due to remaining blends, not included in the 3D-synthetic pro- 
file, NLTE effects that can change the shape of the line profile, 
or the adopted damping constants. Comparison of Figs. Q] and [2] 
shows that the 3D synthetic profiles generally can reproduce the 
observed profiles of the selected clean C i lines somewhat bet- 
ter than the ID synthetic profiles calculated from the HM model 
with a micro-turbulence of l.Okms -1 . 



6.3. The solar carbon abundance 

The final carbon abundance depends only weakly on the assump- 
tion made about Sh- Applying the NLTE corrections to the 3D 
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Fig. 1. Comparison of the observed disc-centre solar spectrum (dashed lines, DI for the two reddest lines and NI for the others) with 
the corresponding 3D synthetic profiles (solid lines) for a selection of clean C i lines. 



LTE results, and computing the average of the abundances of 
Table[3] we obtain: 

A(C)= 8.446 + 0.121 for S H = 

A(C)= 8.498 + 0.110 for S H =1/3 (1) 

A(C) = 8.523 + 0.112 for S H = 1 

For reference, the results for the HM model are: 

A(C)= 8.449 + 0.135 for S H = 

A(C)= 8.503 + 0.116 for S H =1/3 (2) 

A(C) = 8.532 + 0.112 for S H = 1 

The results from CO 5 BOLD and HM models are in very 
good agreement. The carbon abundances from the various lines 
as a function of wavelength are shown in Fig. [3] for the 
CO 5 BOLD model, and in Fig.g]for the HM model. 

Our favoured value is A(C)= 8.498 + 0.110, obtained ap- 
plying the NLTE correction, with Sh = 1/3, to the 3D-LTE 
abundance. If we restrict the abundance determination to the 
lines labelled as "Quality=l" in Table[T| the result is A(C)= 
8.490 ± 0.048. The carbon abundance from the subsample is 
very close to the one obtained from the complete sample, only 
0.008 dex smaller, while the line-to-line scatter is much reduced. 
Nevertheless, we prefer the result from the complete sample, be- 
ing aware that the selection of "good" lines is somehow subjec- 
tive, and hardly changes the mean carbon abundance. 



7. Discussion 

If we consider the EWs of 55 lines in Biemont et al. (1993), 
based on Delbouille disc-centre spectra, together with the log gf 
-values used in that work, we find A(C)=8.518 ±0.137 from our 
3D model. If we use instead the updated log gf -values from 
the NIST database, we obtain A(C)=8.504 ± 0.125. This lat- 
ter value can be com pared to our LTE res ult, based on the 40 
lines in common with Biemont et al. I (fl99l and measured in the 
same solar atlas, of A(C)=8.535 + 0.121 obtained with the 3D 
model, and of A(C)=8.550 + 0.108 obtained with the HM model. 
The LTE abundance based on our complete sample of 98 lines is 
A(C)=8.553 + 0.125 from the 3D model and A(C)=8.569 + 0.1 18 
from t he HM model. Our abu ndance is slightly higher than the 
one of B iemont et al.l (Q993) because of the line selection and 
differences in the EWs for some lines, but the overall agreement 
is very satisfactory. 

3D-NLTE abundances are obtained by applying the 1D- 
NLTE corrections with Sh = 1/3 to the individual 3D-LTE abun- 
dances. All 98 3D-NLTE abundances lie within 3<x of the mean 
value, A(C)= 8.498 + 0.110. Keeping only abundances within 
2cr of the mean value, 94 values meet the cutoff, and the average 
becomes A(C)= 8.492+0.098; within l<x, still 70 values survive, 
and A(C)= 8.485 + 0.049. As expected, the standard deviation 
becomes smaller, but the average remains almost the same. 
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Fig. 2. As Fig.Q] but showing the HM synthetic profile (solid lines) superimposed on the observed solar spectrum (dashed lines). To 
be consistent with Fig.[T| the synthetic profile is computed with Linfor3D, ignoring any blending lines. 



To check the validity of the EW approach, we used line pro- 
file fitting to determine the carbon abundance for two lines which 
are not blended, are weak, and have a very small NLTE cor- 
rection: 538. Onm and 658.7 nm. On average, the result is about 
0.02 dex below the one obtained from the EW. For the majority 
of the selected lines the atomic data of the blending components 
are not very well known. These blending components are, how- 
ever, sufficiently separated that EWs can be measured by fitting 
multiple Gaussian or Voigt profiles. We therefore prefer to use 
EWs, rather than line profile fitting (see Sect.©. 

We find no obvious trend of the abundance with the lower 
level energy, neither for the LTE nor for the NLTE results. 
However, this is not surprising since all carbon lines of our sam- 
ple originate from similar high excitation levels, the range in en- 
ergy being little more than 2 eV. 

There is, however, a clear trend that the 3D-LTE abundances 
increase with equivalent width (see Fig.[5j. The trend in reduced 
or even reversed after application of the NLTE corrections, de- 
pending on the choice of Sh- As illustrated in Fig. [5] there is a 
slight negative trend of the 3D-NLTE abundance with EW for 
Sh = 0.0, and a slightly positive one for Sh = 1/3. The trend 
vanishes for a value of Sh somewhere in the range [0, 1 /3]. The 
corresponding results obtained using the HM model are shown 
in Fig. [6] The behaviour is similar to what is found with the 3D 
model, but the slope of the A(C)-EW relations is systematically 



reduced (more negative) in all cases. The results form the HM 
model suggest that the slope vanished for Sh close to 1 . We note 
that the correlation A(C)lte-EW persists even if we consider 
only weak lines, indicating that the slope is in fact due to NLTE 
effects. 

Without available experimental data on cross sections for 
collisions with neutral hydrogen, one might be tempted to fix 
the value of Sh empirically by requiring a vanishing slope in 
the A(C)-EW plane. Our 3D results would then suggest that 
0.0 < Sh < 1/3. However, we prefer to delay this conclusion 
until a complete 3D-NLTE computation will be available. In the 
meantime, obliged to take a decision, we adopt the intermediate 
value of Sh = 1 /3, which is the favourite value of the Holweger 
school. 

Our carbon abundance s ar e larger than t hose d erived by 
lAllende Prieto et alJ (120021) and lAsplund eTaT] (l2005al) . The dif- 
ference is striking if we consider the [CI] line. While the mea- 
sured EWs are si milar, we adopt an EW which is about 12% 
smaller than that of Asnlund et al. (2005a), due to the correction 
for the blending Fe i line. In spite of this, our d erived A(C) is 
0.1 1 dex higher than that of Asplu nd et alJ d2005al) . This must be 
ascribed to the difference between the CO 5 BOLD solar model 
+ Linfor3D and the hydrod ynamical simu l ation a nd spectral- 
synthesis code employed bv lAsplund et alJ d2005al). The differ - 
ence between the models (see Figure 1 in lCaffau et~ai] d2008)) 
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Fig. 3. The carbon abundance as a function of wavelength of the 
individual C i lines compiled in Table [3] as obtained from the 
3D solar model for different assumptions about the ID NLTE 
corrections, which are all based on the (3D) model. 



affects both the mean temperature structure and the tempera- 
ture fluctuations in the region where the [CI] line is formed. 
Fig. [7] shows that the main contribution to the [CI] line absorp- 
tion comes from the layers between log r = 0, to -2 where the 
differences between two hydrodynamical simulations are largest. 
Si nce the mean te mperature gradient is steeper in the 3D model 
of lAsplund et al.l d2005al) . a lower carbon abundance is needed 
to reproduce the observed equivalent width, in qualitative agree- 
ment with the results mentioned above. The same b ehaviour has 
been noti ced for the [OI] line at 636.3 nm line (see Caffau et al. 
120081 and lAsnlund etai1l2004h . 

The situation is similar if we comp are the average carbon 
abundance. From the permitted C i lines, Asplu ncfet al.l ([2005a) 
derive (A(C)) — 8.36 + 0.03. This value must be compared with 
our value of 8.446 + 0.121, w hich correspond to Sh = 0, con- 
sistent with the assumption of lAsplund et al.l (l2005ah . Part of 
the diffe rence is due to our EW s, which are generally close to 
those of Biem onTet al.l ( 119931) . and therefore larger than those 
of |Xsplund et al. (2005a). Part of the difference is, however, due 
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Fig. 4. As Fig. [3] but showing the carbon abundances obtained 
from the HM model. 



to the different hydrodynamical model and/or the different spec- 
trum synthesis code. 
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Fig. 5. The carbon abundances as a function of EW for the 3D 
model. The two number in the lower part of the plots indicate 
the slope of the best fit linear relation and its 1 cr uncertainty. 
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Fig. 6. The carbon abundances as a function of EW for the HM 
model. 
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Table 3. Line parameters, 3D and ID carbon abundance (LTE) and ID NLTE corrections for our sample of selected C i lines. 



A 


SP 


X 


EW 


log«/ 


Acc. 




A(C) (LTE) 




3D corrections 


1D-NLTE corrections 


nm 




eV 


pm 






3D 


<3D) 


IDlhd 


3D-(3D> 


3D-1D LHD 


1.0 


1/3 


0.0 




872.7126 


DI 


1.26 


0.47 


-8.140 


B 


8.428 


8.401 


8.388 


0.027 


0.040 


0.000 


0.000 


0.000 


872.7126 


NI 


1.26 


0.47 


-8.140 


B 


8.426 


8.399 


8.386 


0.027 


0.040 


0.000 


0.000 


0.000 


872.7126 


KF 


1.26 


0.46 


-8.140 


B 


8.382 


8.360 


8.354 


0.021 


0.028 


0.000 


0.000 


0.000 


872.7126 


NF 


1.26 


0.47 


-8.140 


B 


8.390 


8.369 


8.362 


0.021 


0.028 


0.000 


0.000 


0.000 


477.5907 


DI 


7.49 


1.71 


-2.304 


C 


8.698 


8.709 


8.709 


-0.011 


-0.011 


-0.003 


-0.007 


-0.017 


477.5907 


NI 


7.49 


1.75 


-2.304 


C 


8.711 


8.722 


8.721 


-0.011 


-0.010 


-0.003 


-0.007 


-0.017 


477.5907 


KF 


7.49 


1.42 


-2.304 


C 


8.687 


8.716 


8.697 


-0.029 


-0.010 


-0.006 


-0.013 


-0.030 


477.5907 


NF 


7.49 


1.43 


-2.304 


c 


8.691 


8.720 


8.701 


-0.029 


-0.010 


-0.006 


-0.013 


-0.030 


505.2167 


DI 


7.68 


4.04 


-1.303 


B 


8.428 


8.429 


8.402 


-0.000 


0.026 


-0.007 


-0.015 


-0.036 


505.2167 


NI 


7.68 


4.15 


-1.303 


B 


8.449 


8.449 


8.422 


0.000 


0.028 


-0.007 


-0.015 


-0.036 


505.2167 


KF 


7.68 


3.56 


-1.303 


B 


8.456 


8.482 


8.442 


-0.026 


0.014 


-0.013 


-0.026 


-0.065 


505.2167 


NF 


7.68 


3.54 


-1.303 


B 


8.452 


8.478 


8.438 


-0.026 


0.014 


-0.013 


-0.026 


-0.065 


538.0336 


DI 


7.68 


2.55 


-1.616 


B 


8.444 


8.447 


8.428 


-0.004 


0.016 


-0.007 


-0.014 


-0.034 


538.0336 


NI 


7.68 


2.53 


-1.616 


B 


8.439 


8.442 


8.423 


-0.004 


0.016 


-0.007 


-0.014 


-0.034 


538.0336 


KF 


7.68 


2.19 


-1.616 


B 


8.465 


8.490 


8.458 


-0.025 


0.007 


-0.012 


-0.025 


-0.060 


538.0336 


NF 


7.68 


2.20 


-1.616 


B 


8.468 


8.492 


8.461 


-0.025 


0.007 


-0.012 


-0.025 


-0.060 


658.7608 


DI 


8.54 


1.65 


-1.003 


B 


8.316 


8.322 


8.293 


-0.006 


0.023 


-0.004 


-0.009 


-0.018 


658.7608 


NI 


8.54 


1.68 


-1.003 


B 


8.326 


8.332 


8.303 


-0.006 


0.023 


-0.004 


-0.009 


-0.018 


658.7608 


KF 


8.54 


1.34 


-1.003 


B 


8.343 


8.369 


8.336 


-0.026 


0.007 


-0.007 


-0.016 


-0.030 


658.7608 


NF 


8.54 


1.33 


-1.003 


B 


8.339 


8.365 


8.332 


-0.026 


0.007 


-0.007 


-0.016 


-0.030 


708.7827 


DI 


8.65 


0.72 


-1.442 


C 


8.413 


8.423 


8.401 


-0.010 


0.012 


-0.004 


-0.008 


-0.017 


711.1475 


DI 


8.64 


1.30 


-1.085 


B 


8.348 


8.355 


8.328 


-0.006 


0.021 


-0.004 


-0.009 


-0.020 


711.3180 


DI 


8.65 


2.73 


-0.773 


B 


8.468 


8.467 


8.430 


0.001 


0.039 


-0.005 


-0.011 


-0.026 


711.3180 


NI 


8.65 


2.75 


-0.773 


B 


8.473 


8.472 


8.434 


0.001 


0.039 


-0.005 


-0.011 


-0.026 


711.3180 


KF 


8.65 


2.45 


-0.773 


B 


8.563 


8.582 


8.544 


-0.019 


0.018 


-0.009 


-0.020 


-0.045 


711.3180 


NF 


8.65 


2.39 


-0.773 


B 


8.547 


8.566 


8.529 


-0.019 


0.018 


-0.009 


-0.020 


-0.045 


783.7105 


DI 


8.85 


0.22 


-1.778 


B 


8.363 


8.373 


8.353 


-0.010 


0.010 


-0.003 


-0.006 


-0.012 


801.8564 


DI 


8.85 


0.13 


-2.130 


D 


8.480 


8.489 


8.470 


-0.009 


0.010 


-0.002 


-0.005 


-0.011 


833.5149 


DI 


7.68 


12.30 


-0.437 


B + 


8.606 


8.570 


8.508 


0.037 


0.098 


-0.069 


-0.121 


-0.225 


833.5149 


NI 


7.68 


12.70 


-0.437 


B + 


8.642 


8.605 


8.542 


0.037 


0.100 


-0.069 


-0.121 


-0.225 


833.5149 


KF 


7.68 


10.70 


-0.437 


B + 


8.661 


8.646 


8.602 


0.014 


0.058 


-0.108 


-0.185 


-0.318 


833.5149 


NF 


7.68 


10.60 


-0.437 


B + 


8.650 


8.636 


8.592 


0.014 


0.058 


-0.108 


-0.185 


-0.318 


906.1432 


DI 


7.48 


16.50 


-0.347 


B 


8.643 


8.597 


8.532 


0.045 


0.111 


-0.074 


-0.127 


-0.236 


907.8278 


DI 


7.48 


13.50 


-0.581 


B 


8.640 


8.597 


8.536 


0.043 


0.104 


-0.059 


-0.104 


-0.201 


911.1797 


DI 


7.49 


17.50 


-0.297 


B 


8.665 


8.620 


8.553 


0.045 


0.112 


-0.054 


-0.109 


-0.223 


911.1797 


NI 


7.49 


17.70 


-0.297 


B 


8.678 


8.633 


8.566 


0.045 


0.112 


-0.054 


-0.109 


-0.223 


911.1797 


KF 


7.49 


15.20 


-0.297 


B 


8.715 


8.690 


8.642 


0.025 


0.074 


-0.146 


-0.231 


-0.384 


911.1797 


NF 


7.49 


15.50 


-0.297 


B 


8.739 


8.714 


8.664 


0.025 


0.075 


-0.146 


-0.231 


-0.384 


960.3032 


DI 


7.48 


11.50 


-0.896 


B 


8.727 


8.686 


8.630 


0.041 


0.098 


-0.039 


-0.073 


-0.130 


965.8435 


DI 


7.49 


17.00 


-0.280 


B 


8.559 


8.512 


8.447 


0.047 


0.112 


-0.042 


-0.094 


-0.179 


1012.3871 


DI 


8.54 


11.80 


-0.031 


C+ 


8.682 


8.651 


8.589 


0.031 


0.094 


-0.032 


-0.063 


-0.118 


1012.3871 


NI 


8.54 


11.60 


-0.031 


c+ 


8.665 


8.634 


8.572 


0.031 


0.093 


-0.032 


-0.063 


-0.118 


1012.3871 


KF 


8.54 


9.41 


-0.031 


c+ 


8.677 


8.672 


8.626 


0.006 


0.051 


-0.056 


-0.106 


-0.170 


1012.3871 


NF 


8.54 


9.35 


-0.031 


c+ 


8.671 


8.665 


8.620 


0.006 


0.051 


-0.056 


-0.106 


-0.170 


1054.1241 


DI 


8.54 


2.03 


-1.398 


D 


8.728 


8.725 


8.699 


0.004 


0.030 


-0.009 


-0.017 


-0.042 


1054.1241 


NI 


8.54 


1.95 


-1.398 


D 


8.706 


8.703 


8.677 


0.003 


0.029 


-0.009 


-0.017 


-0.042 


1054.1241 


KF 


8.54 


1.75 


-1.398 


D 


8.800 


8.814 


8.788 


-0.014 


0.011 


-0.016 


-0.032 


-0.074 


1054.1241 


NF 


8.54 


1.65 


-1.398 


D 


8.767 


8.781 


8.756 


-0.015 


0.010 


-0.016 


-0.032 


-0.074 


1068.5345 


DI 


7.48 


20.20 


-0.272 


B 


8.605 


8.556 


8.489 


0.050 


0.116 


-0.067 


-0.117 


-0.205 


1068.5345 


NI 


7.48 


20.00 


-0.272 


B 


8.594 


8.544 


8.478 


0.050 


0.116 


-0.067 


-0.117 


-0.205 


1068.5345 


KF 


7.48 


17.80 


-0.272 


B 


8.668 


8.640 


8.588 


0.028 


0.080 


-0.107 


-0.188 


-0.340 


1068.5345 


NF 


7.48 


17.50 


-0.272 


B 


8.648 


8.620 


8.569 


0.028 


0.080 


-0.107 


-0.188 


-0.340 


1070.7333 


DI 


7.48 


18.50 


-0.411 


B 


8.638 


8.589 


8.524 


0.050 


0.114 


-0.047 


-0.090 


-0.193 


1070.7333 


NI 


7.48 


18.20 


-0.411 


B 


8.619 


8.570 


8.505 


0.050 


0.114 


-0.047 


-0.090 


-0.193 


1070.7333 


KF 


7.48 


16.80 


-0.411 


B 


8.735 


8.707 


8.656 


0.028 


0.078 


-0.127 


-0.207 


-0.357 


1070.7333 


NF 


7.48 


16.00 


-0.411 


B 


8.678 


8.651 


8.601 


0.027 


0.077 


-0.127 


-0.207 


-0.357 


1072.9533 


DI 


7.49 


18.20 


-0.420 


B 


8.644 


8.594 


8.529 


0.050 


0.115 


-0.040 


-0.089 


-0.190 


1072.9533 


NI 


7.49 


17.80 


-0.420 


B 


8.618 


8.568 


8.504 


0.050 


0.115 


-0.040 


-0.089 


-0.190 


1072.9533 


KF 


7.49 


15.80 


-0.420 


B 


8.687 


8.660 


8.611 


0.028 


0.077 


-0.120 


-0.199 


-0.342 


1072.9533 


NF 


7.49 


15.60 


-0.420 


B 


8.672 


8.645 


8.596 


0.027 


0.076 


-0.120 


-0.199 


-0.342 


1075.3985 


DI 


7.49 


4.54 


-1.606 


B 


8.489 


8.467 


8.435 


0.022 


0.055 


-0.012 


-0.025 


-0.055 


1075.3985 


NI 


7.49 


4.48 


-1.606 


B 


8.480 


8.458 


8.426 


0.022 


0.054 


-0.012 


-0.025 


-0.055 


1075.3985 


KF 


7.49 


3.64 


-1.606 


B 


8.477 


8.478 


8.454 


-0.001 


0.023 


-0.023 


-0.045 


-0.095 


1075.3985 


NF 


7.49 


3.61 


-1.606 


B 


8.471 


8.472 


8.449 


-0.001 


0.023 


-0.023 


-0.045 


-0.095 


1174.8220 


DI 


8.64 


14.90 


0.375 


B 


8.368 


8.330 


8.267 


0.038 


0.102 


-0.029 


-0.060 


-0.138 


1174.8220 


NI 


8.64 


14.80 


0.375 


B 


8.362 


8.323 


8.260 


0.038 


0.101 


-0.029 


-0.060 


-0.138 


1174.8220 


NF 


8.64 


12.40 


0.375 


B 


8.392 


8.381 


8.327 


0.011 


0.064 


-0.032 


-0.086 


-0.211 


1177.7540 


DI 


8.64 


6.56 


-0.520 


B 


8.514 


8.488 


8.443 


0.026 


0.071 


-0.015 


-0.030 


-0.073 


1177.7540 


NI 


8.64 


6.58 


-0.520 


B 


8.517 


8.491 


8.446 


0.026 


0.071 


-0.015 


-0.030 


-0.073 


1177.7540 


KF 


8.64 


5.09 


-0.520 


B 


8.491 


8.492 


8.455 


-0.002 


0.036 


-0.026 


-0.053 


-0.125 


1177.7540 


NF 


8.64 


5.09 


-0.520 


B 


8.491 


8.492 


8.455 


-0.002 


0.036 


-0.026 


-0.053 


-0.125 


1184.8730 


DI 


8.64 


5.70 


-0.697 


B 


8.510 


8.490 


8.452 


0.020 


0.058 


-0.012 


-0.025 


-0.060 


1184.8730 


NI 


8.64 


5.75 


-0.697 


B 


8.516 


8.496 


8.458 


0.020 


0.058 


-0.012 


-0.025 


-0.060 


1184.8730 


KF 


8.64 


4.66 


-0.697 


B 


8.534 


8.538 


8.502 


-0.004 


0.032 


-0.023 


-0.049 


-0.108 


1186.2990 


DI 


8.64 


5.81 


-0.710 


B 


8.535 


8.514 


8.476 


0.020 


0.059 


-0.011 


-0.025 


-0.061 


1186.2990 


NI 


8.64 


5.79 


-0.710 


B 


8.532 


8.512 


8.474 


0.020 


0.058 


-0.011 


-0.025 


-0.061 


1186.2990 


KF 


8.64 


4.63 


-0.710 


B 


8.541 


8.544 


8.508 


-0.004 


0.032 


-0.023 


-0.049 


-0.107 


1189.2910 


DI 


8.64 


10.10 


-0.277 


B 


8.518 


8.491 


8.439 


0.027 


0.079 


-0.020 


-0.042 


-0.097 


1189.2910 


NI 


8.64 


10.10 


-0.277 


B 


8.518 


8.491 


8.439 


0.027 


0.079 


-0.020 


-0.042 


-0.097 


1189.2910 


NF 


8.64 


8.45 


-0.277 


B 


8.555 


8.553 


8.506 


0.002 


0.048 


-0.037 


-0.076 


-0.157 


1189.5750 


DI 


8.65 


13.94 


-0.008 


B 


8.530 


8.499 


8.441 


0.031 


0.089 


-0.027 


-0.052 


-0.117 



Table 3. continued. 
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A 


SP 


X 


EW 


log?/ 


Acc. 




A(C) (LTE) 




3D corrections 


1D-NLTE corrections 


nm 




eV 


pm 






3D 


<3D> 


IDlhd 


3D-^3D) 


JU-1ULHD 


1 







1 loV.J /jU 


Nl 






n nnfi 
— u.uufi 


13 


o.jjZ 


8.JU1 


AAl 
S.44Z 


0.031 


0.090 


-0.027 


-0.052 


-0.117 




NF 


8 65 


1 1 45 


—0 008 




8 550 


8 545 


8 49 1 


0.005 


0.058 


-0.023 


-0.073 


-0.169 


1 O^A CM 8ft 


Dl 


8 85 


a en 


—0 565 




8 63^ 


8 604 


8 563 


0.028 


0.069 


-0.008 


-0.017 


-0.045 


1 9S6 91 90 


Dl 


8 8^ 


6 89 


—0 522 




8 597 


8 569 


8 528 


0.028 


0.069 


-0.008 


-0.018 


-0.047 


19S6 0040 


DI 


8.85 


5.75 


—0.598 


[3 


8.535 


8.510 


8.473 


0.025 


0.062 


-0.008 


-0.017 


-0.044 


1258.1590 


Dl 


8.85 


6.00 


-0.536 


B 


8.502 


8.477 


8.439 


0.026 


0.063 


-0.008 


-0.018 


-0.042 


1734.6381 


DI 


9.70 


0.64 


-1.348 


c 


8.352 


8.336 


8.351 


0.016 


0.001 


0.000 


0.000 


0.003 


1744.8600 


DI 


9.00 


21.60 


0.012 


B + 


8.428 


8.393 


8.359 


0.035 


0.069 


-0.006 


-0.015 


-0.058 


1745.5971 


DI 


9.70 


17.40 


0.280 


C 


8.567 


8.534 


8.502 


0.033 


0.066 


0.001 


0.002 


0.008 


1750.5641 


DI 


9.70 


17.50 


0.424 


C 


8.430 


8.396 


8.364 


0.034 


0.066 


0.001 


0.001 


0.008 


1763.7381 


DI 


9.71 


21.00 


0.338 


c 


8.673 


8.639 


8.600 


0.035 


0.073 


0.001 


0.002 


0.013 


1778.9600 


DI 


7.95 


1.40 


-2.246 


B 


8.247 


8.225 


8.239 


0.022 


0.009 


-0.003 


-0.006 


-0.008 


2102.3131 


DI 


9.17 


10.00 


-0.450 b 




8.651 


8.618 


8.580 


0.033 


0.071 


-0.001 


-0.002 


-0.009 


2121.1551 


DI 


9.83 


8.56 


-0.080 6 




8.660 


8.632 


8.599 


0.027 


0.061 


0.000 


0.000 


0.004 


2125.9891 


DI 


9.83 


17.80 


0.490 6 




8.598 


8.566 


8.513 


0.033 


0.085 


0.000 


0.000 


0.004 


3129.7480 


DI 


9.69 


2.60 


-0.570 b 




8.688 


8.673 


8.650 


0.015 


0.039 


0.000 


0.001 


0.002 


2290.6561 


DI 


9.17 


12.15 


-0.182° 




8.793 


8.763 


8.720 


0.030 


0.073 


-0.001 


-0.003 


-0.013 



logg/-values with their quality (column six) are taken from the NIST database; the quality (Acc.) represents the accuracy of the value, ranging from B+, meaning ^j- < 7%, to D, 
meaning ^4- < 50%. log gf with flag b are from Biemont et al. 1 1993), those with flag a are from Asplund et al. 1 2005a) (no quality is indicated in these cases). 
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Table 4. Solar abundance of C in the literature 



A(C) 




Ket 


8.67 


±0.10 


Lambert (1978) 


8.56 


±0.04 


Anders & Grevesse (1989) 


8.58 


±0.13 


Stiirenburs & Holweger (1990) 


8.60 


±0.05 


Grevesse et al. (1991) 


8.60 


±0.10 


Biemont etal. (1993) 


8.55 




Grevesse et al. (1994) 


8.54 




Takeda (1994) 


8.52 


±0.06 


Grevesse & Sauval (1998) 


8.52 


±0.06 


Grevesse et al. (2000) 


8.57 


±0.03 


Holweger (2001) 


8.592 


±0108 


Holweger (2001) 


8.39 


±0.04 


Allende Prieto et al. (2002) 


8.39 


±0.05 


Asplund et al. (2005a) 


8.39 


±0.05 


Scott et al. (2006') 


8.44 


±0.06 


Pinsonneault & Delahave (2006) 


8.39 


±0.05 


Grevesse et al. (2007) 


8.43 


±0.05 


Asplund et al. (2009) 


8.50 


±0.06 


present work 



In Table|4]the carbon abundance determinations in the last 
thirty years are listed. The difference of 0.28 dex from the high- 
est to the lowest value cannot be explained with NLTE effects, 
that, according to our analysis, are about -0.05 dex on average 
for Sh = 1/3. The average of the values in the Table is 8.51 with 
a standard deviation of 0.09. 



8. Conclusions 

Our recommended value for the solar carbon abundance is A(C) 
= 8.50 + 0.06, corresponding to a weak efficiency of the col- 
lisions with neutral hydrogen atoms (Sh = 1 /3), the favourite 
value of the Holweger school. The quoted error is the linear sum 
of a statistical error, 0.02 dex, and a systematic error, 0.04 dex, 
due to the uncertainty of the treatment of the hydrogen colli- 
sions in the NLTE computation. The statistical error has been 
estimated by dividing the line-to-line scatter, 0.11 dex, by the 
square root of the 45 independent lines used in the analysis. 
The value Sh = 1/3 was also ad opted in our invest igations of 
the s olar abundances of o xygen dCaffau et all 12008) and nitro- 
gen (Caffau et all l2009bb . It is not obvious why the same value 
of Sh should apply to different atoms, or even to different lines 
of the same atom. It is comforting that the difference between 
the extreme assumptions about the efficiency of the H collisions 
(Sh = or 1) amounts to only 0.08 dex. 

Our preferred value for the s olar carbon abundance is v ery 
close to the recommendation of iGrevesse & Sauvall (Il998l) . If 
we take this carbon abundance, together with A(N)=7.86 from 
ICaffau et al.l (l2009bh, A(Q)=8.76 from ICaffau et al.1 (l2008h 
and A(Ne)=8.02 (see ICaffau et~alll2009bl for an explanation of 
this choice), we obtain a solar metallicity of Z — 0.0154 and 
Z/X — 0.0 21 1. This value is higher tha n the metallicity recom- 
mended bv lAsplund et al. I (l2005bh and lGrevesse et all (120071) . 
Z = 0.0122, and goes in the direction of reconciling the spectro- 
scopic abundances with the constraints from helioseismology. 

The fact that different 3D hydrodynamical simulations pro- 
vide significantly different results (of the order of 0.1 dex) un- 
derlines the need for further develo pment and validation o f the 
hydrodynamical models. Recently Asplu nd" et al.l d2009l) has 



presented results based on a new generation 3D model, which 
is much closer to the CO 5 BOLD model than the one used by 
Aspl und et al.l d2005al) . and lead to an upward revision of his car- 
bon, oxygen, and iron abundances. It is likely that any residual 
difference between this new result and the present analysis can 
be ascribed to the line selection and to the different treatment of 
hydrogen collisio ns in the NLTE comp utations, although details 
on the analysis o flAsplund et al.l d2009l) are not yet available. The 
excellent agreement of our CO 5 BOLD sola r model with the ob - 
served centre-to-limb variation, shown in Ludwig et al. (2009), 
provides a strong support to the thermal structure of the model 
and to the abundances deduced by its application dCaffau et all 
I2009ah . 
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